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DYNAMIC  STABILIZATION  OF  THE  RAYLEIGH-TAYLOR 
INSTABILITY  ON  LASER-IMPLODED  SHELLS 

I . Introduction 

In  this  paper  a method  for  dynamic  stabilization  of  the  Rayleigh- 
Taylor  instability  on  the  surface  of  an  imploding  fusion  pellet  is 
presented  and  discussed.  The  technique  involves  the  continual  modu- 
lation of  the  driving  laser  beams  so  that  the  potentially  Rayleigh- 
Taylor  unstable  ablator  shell  is  subject  to  a rapidly  and  strongly 
oscillating  acceleration.  A band  of  R-T  modes  can  be  stabilized  by 
this  oscillation  even  though  the  time  average  acceleration  vector 
lies  in  the  destabilizing  direction.  By  adjusting  the  frequency, 
structure,  and  amplitude  of  the  modulation,  the  band  of  dynamically 
stabilized  modes  can  be  made  to  include  the  most  unstable  and 
dangerous  modes.  Thus  considerably  higher  aspect  ratio  shells 
(i.e.,  thinner  shells)  can  implode  successfully  than  had  been  pre- 
viously considered  stable  enough.^-  Both  theory  and  fully  nonlinear 
numerical  simulations'  support  this  conclusion  for  the  case  of 
laser-driven  pellet  implosions.  Similar  beam  modulation  via  charged 
particle  bunching  should  work  to  stabilize  the  most  dangerous 
Rayleigh- Taylor  modes  in  relativistic  electron  beam,  ion  and  heavy 
ion  pellet  implosion  schemes. 

Most  schemes  for  producing  thermonuclear  fusion  by  imploding  a 
solid  sphere  or  hollow  shell  of  D-T  fuel  rely  on  the  symmetry  of  the 
implosion  to  achieve  the  necessary  densities  and  temperatures  at 
peak  compression  to  initiate  the  thermonuclear  burn  process.  The 
prevalent  schemes  suggested  to  date  (laser  beams,  electron  beams, 


Note:  Manuscript  submitted  November  24,  1976. 
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light  or  heavy  ion  beams  incident  on  the  surface  of  a layered  pellet) 
all  involve  deposition  of  energy  in  a thin  ablation  layer  at  the 
surface  of  the  pellet.  This  energy  rapidly  heats  the  dense  shell 
material  causing  a high  pressure  region.  This  pressure  maximum 
drives  the  shell  inward  and  conserves  momentum  by  accelerating  a 


low  density  blow-off  plasma  outward.  Since  the  relatively  low 
density  ablation  layer  is  accelerating  the  higher  density  shell, 
the  system  is  unstable  to  the  Rayleigh- Taylor  modes  which  can  grow 
when  the  acceleration  lies  parallel  to  the  density  gradient. 

In  the  laser-pellet  configuration  the  behavior  of  the  insta- 
bility is  complicated  by  a host  of  important  physical  effects  so 

not  even  the  linear  problem  has  been  solved  adequately.  Computer 

2 3 

codes  such  as  LASNEX  and  FAST2D'  have  been  used  to  study  some  of 
the  properties  of  the  Rayleigh-Taylor  modes  via  numerical  simulation. 
Although  some  of  the  linear  properties  of  the  R-T  modes  have  even 
been  uncovered  analytically,  an  accurate  dispersion  relation 
including  finite  sound  speed,  ablation  layer  gradients,  convection 
and  finite  thermal  conductivity  still  does  not  exist  to  describe 
the  growth  of  these  modes.  Furthermore,  the  nonlinear  growth 
and  turbulence  properties  are  almost  entirely  a mystery. 

Figure  1 shows  the  density,  temperature,  and  pressure  profiles 
at  the  outside  of  a D-T  shell  being  accelerated  by  electron  thermal 
conduction  induced  ablation.  These  profiles  are  presented  exactly 
as  computed  by  the  FAST2D  computer  code  with  the  laser  depositing 
its  energy  at  the  critical  layer  which  lies  off  the  right  hand  side 
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of  the  figure.  The  x = 0 position  was  chosen  to  coincide  with  the 
peak  pressure  point  so  the  blow-off  lies  to  the  right  (x  > 0)  and 
the  accelerating  shell  lies  to  the  left  (x  < 0).  The  ablation  layer 
proper  extends  about  3 or  4 microns  from  the  separation  point  x = 0 
out  into  the  blow-off.  It  results  from  the  strong  nonlinear  depend- 
ence of  the  thermal  conductivity  on  temperature.  A future  note  will 
be  devoted  to  a simple  two-region  quasi-steady-state  analytic  model 
of  these  ablation  layers. 

Our  current  attention  is  directed  toward  studying  the  Rayleigh- 
Taylor  instability  which  occurs  in  the  region  from  -4p  to  Op,, 
indicated  as  Rayleigh-Taylor  unstable,  where  the  density  and  pressure 
gradients  are  opposed.  The  difficulty  in  performing  the  linear 
analysis  arises  from  the  fact  that  the  instability  develops  in  an 
accelerating  frame  of  reference  where  fluid  is  convecting  through 
the  unstable  region  at  several  microns  per  nanosecond. 

II.  The  Numerical  Model 

3 

The  FAST2D  computer  code  solves  the  compressible  ideal  fluid 
equations  in  2D  slab  geometry  with  a strongly  nonlinear  thermal 
conductivity.  Figures  2,  3j  and  4 show  a developing  Rayleigh-Taylor 
mode  in  two  distinct  computer  graphics  displays.  The  (a)  half  of 
the  figures  shows  contours  of  mass  density  at  0.1,  0.2,  0.3,  0.4, 
and  0.5  gm/cc  values.  The  vertical  (y)  direction  in  these  plots 
has  periodic  boundary  conditions  and  the  slab  is  being  accelerated 
to  the  left  in  these  figures  by  the  deposition  of  laser  energy 
at  the  critical  layer  p ~ 0.004  gm/cc  near  the  right  hand  boundary 
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of  the  figure.  In  the  x direction  the  grid  is  stretched  in  the 
low  density  regions  to  the  right  and  left  of  the  slab  proper  to 
move  the  boundaries  and  their  influence  as  far  as  possible  from  the 
R-T  unstable  region  at  the  back  edge  of  the  shell.  The  temperature 
profile  at  the  boundaries  and  at  the  critical  layer  is  programmed  in 
ahead  of  time  and  provides  the  mechanism  by  which  the  various 
modulations  of  acceleration  are  applied.  The  fluid  has  an  outflow 
condition  applied  at  the  right  and  left  boundaries. 

Figures  2b,  3b,  and  4b  show  the  same  data  plotted  perspectively 
with  y increasing  to  the  right  and  x increasing  out  of  the  picture. 
The  density  is  plotted  vertically  at  three  times  during  the  calcu- 
lation. Figure  2 shows  the  initial  conditions,  Figure  3 at  .93 
nanoseconds  shows  the  instability  developing  linearly,  and  Figure  4 
at  1.1  nanoseconds  shows  the  characteristic  nonlinear  development 
of  the  compressible  Rayleigh -Taylor  instability  in  laser-driven 
ablations.  The  nonlinear  properties  of  the  R-T  modes  shown  is  the 
subject  of  future  research.  The  location  x(y)  of  the  p = 0.1  gm/cc 
contour  shown  in  the  figures  is  Fourier  analyzed  to  study  the  spec- 
trum of  the  surface  perturbations  and  to  measure  the  growth  of  the 
Rayleigh- Taylor  unstable  modes.  The  two  main  results  of  the  linear 
Rayleigh- Taylor  analysis  using  FAST2D  are: 

1.  The  long  wavelength  R-T  modes  appear  to  follow  the 
incompressible  y ~ /ka  law  well  beyond  the  wavelengths 
where  sonic  communication  across  a wavelength  takes  longer 
than  a growth  time. 
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2.  There  is  a cutoff  wavelength  at  which  shorter  wavelengths 
are  stabilized  by  the  finite  thickness  of  the  ablation  layer 
and/or  dissipative  terms  and  longer  wavelengths  grow.  The 
dispersion  relation  for  the  modes  will  be  discussed  in 
greater  detail  in  Section  III. 

The  equations  being  solved  in  EAST2D  are  the  following: 


- t^y 


a(pvx) 


S(oV  ) 

y _ 
st 


f— (F  + cV  V ) - |-(dV  V ) 
ax'  x x'  ay'  y x2 * * * * 7 


- f-(oV  V ) - (P  + pV  V ) 
Sx  x y7  dy  y y7 


ar  a n 5/2  sfi 

= _ x £Y  + py  _ r£-'/  2J. 

at  as.  x x ' o axj 


rEV  + PV  - K TP//2  y j 

ay  l y y o ayJ 


-3^ 

where  ~ 10  is  used  for  D-T  and  where  the  temperature  has 

2 2 P 2 

kinetic  units  of  cm  /sec  . That  is,  P = oT  and  E = + IpV  . 

The  value  of  v is  taken  to  be  5/3  throughout  the  material  and 

5/2 

the  T plasma  thermal  conduction  is  used  in  the  shell  as  well. 

Splitting  algorithms  are  used  along  with  the  latest  FCT  algorithms' 
to  give  maximum  speed  on  the  Texas  Instruments  Advanced  Scientific 
Computer  at  '.7RL. 


5 


The  calculations  shown  in  later  sections  were  performed  with 


typical  values  0.5p  * fix  £ l.Oy,,  l.Op,  £ fiy  < 5p  and  with  0.5  pico- 
second < fit  < 2.5  picoseconds.  The  sound  and  flow  speeds  determine 
the  timestep  by  the  usual  Courant  conditions  and  the  sliding  rezone 
capability  of  the  FCT  algorithms  is  used  to  perform  the  calculation 
in  an  accelerating  frame  of  reference  which  keeps  the  dense  shell 
material  centered  in  the  computation  region  throughout  the  run. 

Figure  5 shows  the  temperature  profile  through  the  ablation 
layer  as  calculated  by  FAST2D  at  four  different  times  during  a 
two-dimensional  calculation.  The  solid  line  shows  the  profile 
predicted  using  a simple  two-region  steady  state  theory  which 
inputs  the  numerical  values  for  upstream  density,  temperatures,  and 
burn-through  rate  from  the  FAST2D  calculation.  The  agreement  of 
the  computation  with  the  analysis  is  very  good  considering  the 
rather  coarse  grid  resolution  at  the  "singularity".  Furthermore, 
the  steadiness  of  the  computed  profile  in  time  lends  credence  to 
the  use  of  a steady-state  theoretical  blow-off  model  even  though 
the  ablation  layer  is  only  really  stationary  in  an  accelerating 
frame  of  reference. 

III.  Analysis  of  Dynamic  Stabilization 

The  numerical  results  to  date  are  consistent  with  a very  simple 
dispersion  relation  for  the  growth  of  the  Rayleigh- Taylor  modes, 


2 _ ka  Uk^a 

= !757ir-  3*0  • 


(5) 
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This  has  the  desired  cutoif  at  k^  = k /2}  as  supported  by  simulations, 
and  has  V ~ */£a  for  long  wavelength,  also  supported  by 

simulation.  In  Fq.  (5)  the  acceleration  of  the  shell  is  a,  k is  the 

P 

inverse  scale  length  of  the  ablation  layer,  and  k = ^ is  the  wave 

A. 

number  of  the  Rayleigh-Taylor  modes  in  the  y-direction  along  the 

shell.  Defining  f = k/k  gives 

0 


2 


V 


k a 
P 


f e 
: 1+e 


(6) 


whc  the  shortest  linearly  unstable  mode  has  X = 4n/k  . Figure  6 

P 

depicts  this  phenomenological  dispersion  relation,  plotting  y/(k  a) 

P 

versus  ohe  nondimension  wave  number  e over  the  linearly  unstable 
range  0 £ e £0.5.  Maximum  growth  occurs  at  e ~ £ (more  correctly 
~ 0.243)  and  the  value  of  the  growth  rate  maximum  is  roughly 


max 


(7) 


The  first  term  in  Eq.  (5)  and  (6)  encompasses  the  usual  simple 
Rayleigh-Taylor  growth  without  infinite  growth  at  short  wavelength. 
The  second  term  takes  the  growth  to  zero  but  is  reactive  in  this 
simple  form--not  dissipative. 

In  the  run  shown  in  Figure  1 and  in  Figure  5,  an  estimate  of 

k can  be  made.  Taking  1.5p  = l/k  as  the  scale  length  for  the  high 
P P 

density  half  of  the  ablation  layer,  the  cutoff  wavelength  from  the 
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simple  dispersion  relation  is 


X = 2tt  x 2/k  = 6n p, 
c ' p 


(8) 


in  good  agreement  with  the  observation  that  mode  4 (X^  = 25u) 
grew  rapidly  but  mode  5 ( X^  = 20p)  did  not. 

To  study  the  effects  of  modulating  the  laser  beam  numerically, 
I varied  the  plasma  temperature  at  the  critical  surface  according 
to  a prescribed  functional  form,  constant  plus  a sinusoidal 
oscillation.  To  study  the  same  phenomenon  analytically  in  an 
illuminating  but  approximate  manner,  the  dispersion  relation 
Eq.  (5)  was  converted  to  a second-order  ordinary  differential 
equation  for  the  perturbation  mode  amplitude  in  which  the  accele- 
ration term  is  varied  parametrically.  This  equation, 


sf*  . r y(*> 
at2  L 1+e 


e 


4e2k 

Y<a(t))]A, 


(9) 


like  the  dispersion  relation  on  which  it  is  based,  has  only  a 
phenomenological/intuitive  basis  because  the  linear,  time-constant- 
coefficient  problem  hasn't  been  solved  satisfactorily  to  date.  The 
acceleration  in  the  cutoff/ dissipation  term  is  taken  as  a time 
average  although  some  additional  stabilization  might  result  from 
including  in  the  theory  the  oscillating  thickness  of  the  ablation 
layer  itself. 
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Let  a(t)  = aQ  + a^  cos  urt;  where  aQ  is  the  time  average  part  and 
a^  is  the  amplitude  of  the  oscillating  acceleration.  The  ordinary 
differential  equation  (9)  becomes 


— ~ = r r—~  (k  a + k a,  cos  ujt) 
2 L 1+e  o o pi 


where  m is  the  frequency  of  the  modulation.  The  time  dependence 
is  assumed  here  to  be  simple,  periodic,  and  sinusoidal.  Later  I will 
discuss  pulse  shape;  there  is  no  reason  to  assume  that  the  sinus- 
oidal modulation  is  optimal  but  at  least  it  is  analytically  tractable. 
Identifying  with  the  canonical  Mathieu  Equation  as  treated  in 
Morse  and  Feshbach  (see  Figure  7)  requires  20  = tut  and  gives 

2 

.2.  ..  4k  an  e 4k  ea  16k  e 

H - r =°»  2« + -y-*  — §-  ao  1 *•  (u> 

50  it  {1+0)  u;  (1+e)  3w  J 

We  identify  coefficients  in  the  canonical  Mathieu  Equation, 


and  b 


^"6,1  - -y-(l+e)^ 


(14) 


are  obtained  as  the  appropriate  non-dimensional  Mathieu  Equation 

coefficients.  Here  both  a^  and  a^  are  positive  constants. 

With  the  values  of  h and  b determinable  from  Eqs.  (l4) 

above>we  turn  to  Fig.  7 which  shows  the  stability  regions  for 

solutions  of  Mathieu' s Equation.  The  regions  of  stability  and 

instability  are  separated  by  the  solid  lines  in  the  figure  end  we 

2 

need  consider  only  the  upper  half  plane  because  h is  determined 

from  the  physical  parameters.  Above  and  to  the  left  of  the  line 

labelled  b at  least  one  of  the  two  independent  solutions  of 
L 

Eq.  (ll)  grows  without  bound.  Two  stable  bands  are  shown,  one 
between  the  lines  b^  and  bg  and  the  other  in  the  lower  right  hand 
corner  for  large  positive  b.  The  first  region  is  accessible  by 
oscillating  the  laser  intensity  but  the  second  is  not  because 
b is  generally  less  than  h /2.  This  limit  is  shown  as  the  dashed 
line  in  Fig.  7.  For  finite  values  of  6 the  locus  of  (b,h)  values 
as  a function  of  e from  Eq.  (l4)  is  displaced  to  the  left  from 
the  dashed  line  at  sufficiently  short  wavelengths.  Under  such  a 
displacement  only  a band  of  Rayleigh- Taylor  modes  is  stabilized 
and  we  need  to  study  the  extent  of  this  stable  band. 

The  first  stability  boundary  is  given  by 


b 


L 


+ 7h8 

32  1024x32 


(15) 


10 


i 


J 


as  determined  from  Abramowitz  and  Stegun  . From  the  same  source 
the  second  stability  boundary  is  found  to  be 


The  subscripts  L and  S refer  to  the  long  and  short  wavelength 
limits  of  the  first  stable  region  respectively.  We  want  to  find 
the  broadest  band  of  stabilization  that  can  be  achieved  given  some 
control  over  uj  and  a^  for  the  modulation. 

One  needs  to  determine  first  the  value  of  f for  which  the 
linearly  most  unstable  mode  e ~ £ lies  right  on  the  short  wave- 
length stability  boundary.  We  do  this  by  setting  b(e  = h)  = b (h) 

O 

2 

and  solving  for  the  value  of  h . The  result  is  relatively  simple 
for  small  6(i.e.,  aQ/a^  < u) • 

hs  “ 1v/<1-Tr>2  * ’ - <17> 

Using  h = 2 in  Eq.  (l4)  allows  us  to  determine 

O 

4k  a. 

f1  = — = 10  (18) 

uu 

as  the  required  combination  of  k^,  a^  and  uo  such  that  the  linearly 
most  unstable  mode  is  just  stabilized  by  the  modulation.  Perhaps 
f can  be  increased  beyond  this  limit  even  though  a narrow  band  of 
short  wavelength  unstable  modes  might  appear  but  f = 10  is  a 
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simple,  easy  to  remember  criterion. 

Given  6 and  f^,  we  can  solve  for  the  long  wavelength  limit 
of  the  dynamically  stabilized  band  by  setting  b = bT  and  solving 
for  the  value  at  which  marginal  dynamic  stability  occurs. 

We  get 


which  has  the  approximate  solution 


The  width  of  the  stabilized  band  is  therefore  given  by  the  simple 
formula 


(21) 


Clearly  the  larger  the  ratio  of  a^  to  aQ,  the  longer  the  wavelengths 
that  can  be  dynamically  stabilized.  In  an  earlier  computer  calcu- 
lation RUN  1G  (see  next  section  and  Figure  8 ) the  ratio  6 = a^a^ 

was  measured  to  be  of  the  order  \ to  It  is  not  clear  what  to 

? 10 
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use  as  a^  because  of  the  strong  spatial  gradients  in  the  region 


prone  to  Rayleigh-Taylor  instability  . Thus  k k as  an 

L y o 


estimate.  Since  the  ablation  layer  thickness  k lay  in  the  range 

P 


lp  to  2p  (also  estimated  by  seeking  the  most  unstable  wavelength 


X ~ 13p  and  the  cutoff  wavelength  X < 10p),  the  longest  stable 
m c 


wavelength  for  that  run  could  be  expected  to  lie  in  the  ballpark 


range  of  30p  to  70  p.  In  fact  the  mode  with  X = 60p  grew  slowly. 


To  maximize  the  effectiveness  of  this  dynamic  stabilization 


one  clearly  wants  to  minimize  £ in  Eq.  (20).  This  is  accomplished 

L 


by  reducing  6 and  increasing  f^  subject  to  the  constraint  that 


short  wavelengths  "uncovered"  at  large  k by  pulling  the  dynamic 


stabilization  "blanket"  to  small  k (minimizing  £ ) do  not  fracture 

1j 


the  shell.  There  is  also  another  constraint  on  increasing  f^  which 


arises  because  and  a are  really  coupled.  When  uu  is  too  small 


(to  maximize  f^),  the  value  of  a^  also  decreases  because  the  ablation 


layer  and  the  Rayleigh-Taylor  unstable  region  beneath  (see  Figure  l) 


can  adjust  rather  quickly  to  changes  in  the  driving  pressure. 


Assuming  that  the  Rayleigh-Taylor  unstable  region  is  oscillated  in 


position  according  to  l ~ + 2^  cos  uot  , the  acceleration 


a^  ~ £^uu  and  is  proportional  to  frequency  squared.  As  uu  decreases, 


at  least  for  small  uu,  so  does  a.. . Thus  fn  « 4k  jL  is  roughly 

’ 1 1 p 1 


independent  of  uu,  if  is  held  constant  as  uu  is  varied. 


Further  work  is  definitely  required  to  pin  down  more  accurately 


what  the  accessible  regime  of  values  of  f is.  We  turn  now  to  the 
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very  important  question  of  how  this  dynamic  stabilization  effect 
scales  with  size.  Larger  pellets  with  thicker  shells  implode  more 
slowly  so  it  is  correspondingly  important  to  stabilize  proportion- 
ately longer  wavelengths.  When  the  ablating  shell  thickness  AR  is 
increased,  aQ  decreases  as  AR  1 since  the  average  shell  acceleration 
goes  inversely  as  the  total  mass  at  a fixed  pressure.  The  oscillating 
acceleration  a^  is  unaffected  by  the  shell  thickness  so  the  ratio 
6 = aQ/a^  decreases  as  AR  \ This  has  two  advantages  when 

scaling  up  to  larger  pellets.  First,  the  long  wavelength  stable-to- 
unstable  crossover  point  at  kT  « 86k  /f  moves  progressively  toward 

L pi 

longer  wavelength  (smaller  e)  as  6 approaches  zero.  This  means 
that  longer  wavelengths  can  be  stabilized  by  a given  modulation 
frequency  and  amplitude  in  large  pellets  than  in  small  pellets. 

Second,  a smaller  value  of  6 means  that  correspondingly  lower 
modulation  frequencies  can  be  tolerated  since  a^  > aQ  is  easier 
to  ensure  in  large  pellets  than  in  small. 

I 

The  conclusion  is  that,  by  varying  uo,  the  most  unstable  mode 

with  wave  number  k = £k  can  be  moved  right  to  the  boundary  of 

m p 

the  stable  region  near  b ~ 1.8,  h ~ 1.9  (or  beyond)  to  maximize 
the  extent  of  the  dynamically  stabilized  portion  of  the  Rayleigh- 
Taylor  spectrum.  Furthermore,  as  will  be  seen  through  complete 
numerical  simulations  in  the  next  section,  the  region  actually 
stabilized  is  quite  broad  enough  to  increase  the  workable  aspect 
ratio  of  the  imploding  shell  considerably. 


Ih 


IV.  Numerical  Simulations  of  Dynamic  Stabilizations 

Several  sets  of  computer  runs  have  been  performed  to  test 
dynamic  stabilization  in  a realistic  system  which  contains  the 
finite  sound  speed/shock  properties  of  the  fluid  involved  and  which 

also  represents  accurately  in  two  dimensions  the  nonlinear  thermal 

•5 

conductivity  effects  expected.  The  basic  code,  IAST2D-,  has  been 
described  elsewhere.  The  calculations  are  performed  on  a two- 
dimensional  Cartesian  domain  with  periodic  boundary  conditions  in 
the  y direction.  The  laser  is  assumed  to  be  incident  from  the 
x = boundary  of  the  system  and  deposits  energy  at  the  critical 

density  p ~ 0.00U  gm/cc  in  such  a way  that  the  temperature  at  the 
critical  density  is  specified  as  a controlled  time  dependent  para- 
meter. The  location  of  the  critical  layer  is  allowed  to  change 
naturally  with  time  and  the  translation  velocity  of  the  coordinate 
system  is  chosen  so  that  the  dense  shell  material  remains  in  the 
center  of  the  region  of  computation. 

Figure  8 summarizes  the  results  of  a set  of  early  calculations 
performed  on  a 40  x 40  finite  difference  grid.  Three  different 
runs  were  performed  with  differing  laser  intensity  profiles.  The 

reference  run,  labelled  RUN  1C,  began  with  a temperature  at  the 

lU  2 2 

critical  layer  of  1 x 10  cm  /sec  . This  corresponds  to  about 

l4  2 2 

175  ev.  The  temperature  linearly  increased  to  2 x 10  cm  /sec  in 
0.5  nanosecond  and  then  was  held  constant  at  that  level  for  the 
remainder  of  the  run.  The  solid  curve  in  the  figure  shows  the 
velocity  of  the  shell  as  a ^unction  of  time  up  to  the  point  at 
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which  the  shell  buckled  due  to  the  Rayleigh-Taylor  instability. 

There  was  no  laser  modulation  in  the  reference  calculation  and  the 

time  of  buckling  was  determined  to  be  that  time  when  the  maximum 

perturbation  of  the  ablation  layer  was  ± 2p,.  This,  of  course, 

corresponds  to  a far  smaller  perturbation  of  the  inside  edge  of  the 

shell.  Figure  3 shows  the  buckling  point  pictorially. 

The  calculation  labelled  IB  in  the  figure  (x's)  was  performed 

using  the  reference  temperature  profile  plus  a 25$  oscillation  of 

the  laser  intensity  superimposed.  The  amplitude  of  the  oscillation 

2 2 

therefore  was  5 X 10  cm  /sec  and  the  oscillation  frequency  was 
taken  to  be  three  cycles  per  nanosecond.  The  shell  in  RUN  IB  accele- 
rated somewhat  faster  than  the  reference  case  1C  because  the  extra 
energy  deposited  at  the  peak  of  the  oscillation  exceeds  the  energy 
decrement  in  the  trough  of  the  oscillation  because  the  thermal 
conductivity  is  nonlinear.  The  shell  also  accelerates  for  a some- 
what longer  period  of  time,  reaching  a terminal  velocity  roughly 
twice  the  velocity  attained  in  the  reference  case.  The  calculation 
was  terminated,  not  when  the  shell  buckled,  but  when  the  laser 
pulse  burned  through  the  ablator  material.  Thus  the  shell  did  not 
really  go  Rayleigh-Taylor  unstable  at  the  run's  termination  but 
rather  was  simply  too  thin. 

RUN  1G,  plotted  as  triangles  in  Figure  8,  was  performed  with  a 

50$  oscillation  at  2.25  cycles/nsec  superimposed  on  the  linear  ramp 

plus  constant  temperature  profile.  The  oscillating  tempera- 

l4  2 2 

ture  at  the  critical  layer  was  ± 10  cm  /sec  and  the  acceleration 
is  correspondingly  faster  than  either  the  reference  case  1C  or  the 
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2%  oscillation  case  IP.  The  shell  again  burns  through  rather  than 
buckling.  The  velocity  attained  before  bumthrough  is  now  roughly 
2.5  X 10  cm/sec  and  more  than  six  times  the  kinetic  energy  per 
particle  has  been  translated  to  the  shell  in  the  modulated  case  as 
opposed  to  the  unmodulated  reference  case.  The  time  of  acceleration 
and  burnthrough  for  the  5 Of  modulation  in  RUM  1G  is  somewhat  shorter 
than  the  25$  oscillation  because  the  rate  of  acceleration  and 
ablation  of  material  is  somewhat  larger. 

In  the  reference  case  1C  the  linear  growth  rates  for  the  modes 
near  the  peak  of  Figure  6 could  be  determined  roughly  from  the 
calculation  and  supported  at  least  semi-quantitatively  the  dispersion 
relation  of  Eq.  (5).  Detailed  linear  numerics  is  very  difficult  to 
perform  and  analyze  because  the  "equilibrium"  is  constantly  changin'- 
in  time  and  because  the  propagation  of  large  amplitude  compressional 
and  sound  waves  through  the  shell  clearly  affects  the  Rayleigh- 
Taylor  growth  strongly.  Sound  waves,  for  example,  can  locally 
reverse  the  density  gradients  and  temporarily  stabilize  the 
Rayleigh -Taylor  modes  over  half  a cycle.  Of  course,  the  modes  are 
presumably  more  unstable  in  the  second-half  cycle  of  a sound  wave 
when  the  gradients  become  steeper  but  there  is  no  reason  to  expect 
the  phenomena  to  cancel  on  average.  The  finite  amplitude  sound 
waves  also  effectively  transmit  perturbations  at  the  outer  edge  of 
the  shell  throughout  the  shell  volume  during  the  initial  ramp-up 
phase . 

To  pursue  dynamic  stabilization  further,  runs  with  a thicker 
shell  were  clearly  required 
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and  runs  with  by  > bx  were  needed  to  study  the  effects  of  dynamic 

stabilization  on  longer  wavelengths.  Greater  resolution  was  also 

called  for  to  increase  the  breadth  of  the  potentially  unstable 

portion  of  the  spectrum.  A set  of  two  calculations  called  DUMP2A 

were  performed  on  an  80  x 80  finite  difference  grid  with  a somewhat 

higher  base  temperature  and  with  a shell  thickness  of  37p  initially. 

The  idea  was  to  perform  an  entire  reference  calculation  without 

modulation  and  then  to  restart  that  identical  calculation  in  the 

middle  when  initial  wave  amplitudes  had  grown  out  of  the  early 

random  noise  and  mode  growth  was  beginning  to  occur  in  an  organized 

linear  manner.  The  modulated  case  was  restarted  with  a 33 

oscillation  at  two  cycles/nsec  and  the  evolution  of  the  same 

initial  conditions  with  and  without  modulation  is  compared.  The 

l4  2 2 

driving  temperature  ramped  linearly  up  to  3 X 10  cm  /sec ' and 
then  was  held  constant  for  the  duration  of  the  reference  calculation. 
The  amplitude  of  the  oscillation  added  to  this  constant  temperature 
was  1 x 10  cm  /sec  . Earlier  Figures  1 and  5 showed  the  ablation 
layer  details  and  theoretical  computations  for  the  ablation  layer 
structure  for  the  DUMP2A  run.  The  data  in  Figure  9 are 

plotted  for  Modes  1,  2,  and  4 with  wavelengths  of  100^,  50^,,  and 
25p.,  respectively.  From  about  3*25  nsec, on  roughly  exponential 
growth  can  be  seen  for  these  three  modes  in  the  reference  case. 

The  problem  was  initialized  with  modes  4 and  8 excited. 

At  3*25  nsec  the  mode  structure  was  sufficiently  well  developed 
and  yet  the  amplitudes  were  sufficiently  small  that  the  modulation 
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could  be  started  and  the  dynamic  stabilization  properties  of  this 
modulation  could  be  tested  in  a difficult  case  which  had  been 
allowed  to  grow  linearly  for  more  than  3 nsec.  The  two  data  curves 
shown  in  Figure  9 apply  to  modes  2 and  4 and  demonstrate  that  the 
amplitude  of  the  ablation  layer  perturbations  became  essentially 
oscillatory  and  stopped  growing  at  a couple  tenths  of  a micron 
amplitude.  In  the  reference  calculation  the  shell  became  Rayleigh- 
Taylor  unstable  and  broke  at  about  6-6.5  nsec.  The  two  micron 
perturbation  limit  was  again  used  as  a signature  of  shell  breaking 
for  the  unmodulated  reference  case  and  the  modulated  case  showed 
shell  burnthrcugh  rather  than  Rayleigh-Taylor  buckling  (Fig.  3,4). 

It  is  important  to  note  that  often  twines  the  longest  wavelength 
in  the  system  acquires  appreciable  energy  because  the  system  computed 
does  not  provide  a continuum  of  modes  near  the  longest  mode.  In 
DUMP2A,  for  example,  the  unmodulated  case  showed  mode  1 almost  as 
large  as  mode  2 and  mode  4 even  though  the  linear  growth  of  mode  1 
should  be  the  slowest  of  the  three  (or  at  least  slower  than  mode  2). 
The  only  solution  to  this  problem  of  k-space  discreteness/blocking 
is  to  use  finer  grids  and  larger  systems. 

Figures  10  and  11  summarize  results  from  two  runs  performed 
to  study  the  long  wavelength  effects  in  a system  with  the  thicker 
37  p shell.  These  cases,  RUN  2B  and  RUN  2C,  used  the  same  reference 
temperature  profile  as  DUMP2A  but  the  temperature  oscillation 
amplitude  was  increased  to  50 1o  as  shown  in  Figure  10a.  The  period 
of  the  oscillation  was  reduced  to  1 cycle/nanosecond  to  maximize  f^ 


and  hence  to  broaden  the  band  of  dynamic  stabilization  as  much  as 
possible  commensurate  with  keeping  the  shell  entirely  within  the 
region  of  computation.  Figure  11a  summarizes  the  average  shell 
velocity  history  for  these  runs  which  differed  only  in  the  value  of 
6y  (6y  = 1.25m-  for  RUN  2E  and  6y  = 3.75p  for  RUN  2C). 

Again  the  shell  burned  through  rather  than  breaking  due  to 
Rayleigh- Taylor  instability.  iTie  shell  attained  a velocity  of 
3.5  x 10  cm/sec  and  showed  no  serious  buckling  before  burnthrough. 
The  outer  surface  perturbations  for  the  two  runs  for  the  three 
longest  wavelength  modes  are  plotted  in  Figures  10b  and  lib. 

In  both  cases  the  dashed  line  shows  the  linear  growth  pattern 
expected  according  to  Eq.  (5)  for  the  longest  wavelength  in  the 
system  (slowest  growth).  In  both  RUN  2B  and  in  RUN  2C  the  mode 
growth  falls  far  below  the  linear  growth  lines  even  though  two 
out  of  the  three  modes  plotted  should  grow  much  faster.  In  both 
cases  the  outer  surface  perturbation  at  burnthrough  was  at  most  a 
couple  tenths  of  a micron,  a good  order  of  magnitude  below  the 
level  used  to  signal  Rayleigh- Taylor  buckling  in  reference 
calculations  such  as  shown  in  Figure  3. 

Figure  lib  shows  that  modes  of  at  least  300g.  wavelength  along 
the  shell  are  dynamically  stabilized  by  modulating  the  laser 

7 

intensity.  Using  a time  of  6nsec  and  a speed  of  3 X 10  cm/sec 

15  2 

gives  an  average  acceleration  of  aQ  ~ 5 X 10  cm/sec  for  the  run. 
The  fluctuating  acceleration  also  shows  in  Figure  11a  and  is  2-4 
times  the  average.  Thus  a ~ 1 - 2 x 10^  cm/sec^.  Plugging  into 
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Eq.  (l8)  defining  f^  and  using  k^1  ~ 1.5p  and  it  = 2tt  x 10^  sec  1 
gives  f ~ 7 - 13  which  nicely  brackets  the  theoretical  target 
value  of  f^  ~ 10. 

In  6 nsec  at  an  average  velocity  of  roughly  1.5  x 10  cm/sec 
the  shell  travels  900^.  If  the  initial  spherical  shell  were  of 
order  1000 p.  in  radius,  this  case  would  represent  a (relatively) 
successful  implosion  of  a 30  : 1 aspect  ratio  system. 

V . Conclusions 

The  theory  and  calculations  presented  in  this  paper  support  the 
conclusion  that  an  appreciable  band  of  the  unstable  Rayleigh- Taylor 
spectrum  on  the  surface  of  an  imploding  laser  pellet  can  be  stabi- 
lized by  modulating  the  intensity  of  the  driving  laser  beams.  By 
just  stabilizing  the  linearly  most  unstable  mode,  the  band  of  complete 
stabilization  can  be  extended  to  modes  with  wavelengths  more  than 
ten  times  the  shell  thickness.  'The  growth  of  slower  long  wavelength 
modes  is  not  stabilized  by  the  modulation  but  the  implosion  time  is 
sufficiently  short  that  tc.ese  residual  modes  should  not  have  time  to 
seriously  distort  the  shell  before  implosion  occurs. 

The  crucial  question  of  how  large  the  constant  f in  l.q.  (lS) 
can  be  made  and  hence  how  broad  a band  of  modes  can  be  stabilized  is 
at  best  partially  answered.  In  one  of  the  early  Uo  x 40  calculations 
a frequency  of  6 cycles/nsec  was  tried  and  the  shell  went  Rayleigh- 
Taylor  unstable  almost  as  if  no  modulation  had  been  applied.  A high 
frequency  cutoff  like  this  is  expected  because  temperature  oscillations 
at  the  critical  layer  can't  be  felt  at  the  ablation  layer  unless 
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their  period  is  longer  than  the  heat  diffusion  time. 

A kO  x UO  run  was  also  performed  with  a modulation  frequency  of 
1 cycle/nsec  but  this  proved  to  be  too  long  on  such  a small  system. 

In  the  troughs  of  the  oscillation  the  shell  tried  to  expand  out  of 
the  computation  region,  a problem  later  rectified  in  a series  of 
90  x 80  calculations.  For  the  typical  parameters  of  interest  here 
frequencies  of  ~ 0.5  cycle/nsec  and  slower  have  yet  to  be  considered 
and  may  work  best  of  all.  Clearly,  however,  the  period  of  the 
modulation  can't  be  longer  than  a few  growth  times.  It  may  even  be 
possible  to  leave  short  wavelength  modes  unstable  to  stabilize  even 
longer  modes.  Certainly  the  range  of  interesting  oo  values 
spans  about  an  order  of  magnitude. 

Dynamic  stabilization  works  in  the  laser-driven  implosion  because 
the  shell  must  compress  somewhat  to  transmit  the  acceleration  from 
the  ablation  layer  on  the  outside  throughout  the  shell.  Therefore 
an  outward  acceleration  is  necessary  to  relieve  the  shell  compression 
when  the  laser  intensity  reaches  a minimum.  The  fluid  compressibility 
acts  as  a spring  so  that  the  fluid  acceleration  seen  in  the  Rayleigh- 
Taylor  unstable  region  actually  changes  sign  even  though  the  laser 
beam  pushes  the  shell  only  inward.  Tor  dynamic  stabilization  to 
work,  an  average  destabilizing  acceleration  must  at  least  change 
sign  and  the  oscillating  component  must  exceed  appreciably  the 
constant  component  in  amplitude. 

The  theory  was  performed  using  the  Mathieu  Equation  which  in 
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turn  assumes  a sinusoidal  variation  of  the  acceleration.  Work  is 


i 


currently  underway  to  provide  a numerical  capability  for  integrating 
the  time  dependent  mode  amplitude  equations  for  a spectrum  of 
Rayleigh- Taylor  modes  from  k = to  k = O.Olk^.  Then  an  approxi- 
mate evaluation  of  the  relative  efficiency  of  various  non-sinusoidal 

l 

acceleration  profiles  can  be  carried  out  rather  quickly.  A good  ID 
piggyback  code  would  have  to  be  used  to  study  in  detail  the  exact 
amount  of  stabilization  provided  by  laser  modulation  with  any  given 
profile . 

Since  the  portion  of  the  laser  energy  which  gets  deposited  at 
the  critical  surface  depends  on  a number  of  laser  and  fluid  para- 
meters and  since  the  thermal  conduction,  which  converts  temperature 
at  the  critical  layer  to  acceleration  at  the  ablation  layer,  is 
highly  nonlinear,  one  can  only'  guess  what  the  laser  intensity  pro- 
file would  have  to  be  in  order  for  the  acceleration  to  be  sinusoidal 
as  assumed  in  the  theory.  It  is  safe  to  guess  that  the  numerical 
simulations  performed  to  date  are  not  yet  optimum  and,  therefore, 
the  stabilization  conclusions  presented  are  correspondingly  con- 
servative. What  also  remains  to  be  determined  is  the  effect  on 
dynamic  stabilization  of  the  host  of  non-ideal  effects  such  as 
radiation,  atomic  physics,  realistic  equations  of  state,  self- 
magnetic  fields,  and  supra- thermal  particle  effects  which  can  also 
broaden  the  ablation  layer  and  hence  change  the  short  wavelength 
cutoff. 

In  pellet  fusion  schemes  where  the  implosion  requires  electron, 
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ion  or  heavy  ion  beams  to  provide  the  energy  required  at  the 
pellet  surface,  the  ablation  layer, and  potentially  R-T  unstable 
region  beneath,  nevertheless  acts  as  a spring  storing  up  compressional 
energy.  If  the  impinging  particle  beam  were  interrupted  or  reduced 
in  intensity,  the  stored  compressional  energy  is  again  available  to 
provide  a reverse  acceleration  just  as  in  the  case  of  a modulated 
laser.  Thus  dynamic  stabilization  is  potentially  useful  for  these 

charged  particle  beam  fusion  schemes  as  well.  A shaped  cavity 

5 

could  be  used  to  modulate  relativistic  electron  beams  and  similar 
electromagnetic  or  magnetic  techniques  during  the  formation,  accele- 
ration and  transport  of  ion  beams  might  readily  be  employed.  Tn  all 
these  cases,  based  on  our  laser-DT  shell  calculations,  modulation 
frequencies  of  a few  tenths  to  a few  cycles  per  nanosecond  seem  to 
be  required.  This  would  demand  particle  bunches  from  a few  inches  to  a 
few  feet  in  length.  For  large  charged-particle -driven  systems  one 
might  expect  modulations  even  slower  than  one  cycle  per  nanosecond 
to  be  effective  because  the  sonic  response  time  of  a thicker  ablation 
layer  would  be  correspondingly  slower.  It  is  this  ablation  layer  response 
which  determines  the  oscillating  component  of  the  acceleration  and 
hence  the  width  of  the  dynamically  stabilized  band  of  Rayleigh- Taylor 
modes. 
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0.5  gm/cc  are  plotted  here  and  in  subsequent  figures. 
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to  signal  shell  failure,  about  2^  in  the  largest  harmonic 
on  the  outside  of  the  shell.  The  perturbation  of  the 
inner  surface  is  much  less. 


cm 


X 

Figure  5 


Comparison  of  ablation  layer  structure  computed  by 
FAST2D  (data  points)  with  a simple  two-region  analytic 
model  (solid  curve).  The  analytic  model  solves  the 
three  time  independent  conservation  equations  to  the 
right  of  the  inflection  point  and  the  time  independent 
constant  density  thermal  conduction  problem  to  the  left. 


Figure  6 

Phenomenological  dispersion  relation  for  ablation 
Rayleigh- Taylor  modes,  kp  is  l/L^  where  Lp  is  the 

characteristic  length  in  the  ablation  layer,  a is  the 
acceleration  and  y the  mode  growth  rate.  Long 
wavelength  growth  rates  are  /ka  and  the  short  wave- 
lengths stabilize  due  to  the  cutoff  at  k = kp/2. 

Maximum  growth  occurs  at  k « k /U  and  is  about  one-third 
of  the  simple  a value.  p 
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Figure  7 


Stability  diagram  for  the  Mathieu  Equation.  The  region 
of  b-h  parameter  space  between  the  curves  labelled  b^ 
and  b is  dynamically  stabilized  by  the  oscillatory 

modulation.  The  dashed  curve  is  the  locus  of  the 
Rayleigh-Taylor  spectrum  in  the  limit  of  large  amplitude 
oscillations.  The  locus  displaces  to  the  left  slightly 
as  the  average  acceleration  approaches  the  size  of  the 
oscillating  component. 
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Figure  8 

Comparison  of  shell  velocity  versus  time  trajectories 
for  three  calculations  demonstrating  dynamic  stabilization 
of  the  laser  shell  Rayleigh- Taylor  mode.  With  dynamic 
stabilization  the  shell  does  not  buckle  due  to  the 
Rayleigh- Taylor  instability  before  the  ablator  material 
simply  burns  through. 
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Figure  10 

Dynamic  stabilization  of  intermediate  wavelength  Rayleigh 
Taylor  modes  for  a 1 cycle  per  nanosecond  modulation. 

The  upper  curve  gives  the  critical  layer  temperature 
history  for  the  calculation.  The  data  in  the  lower 
figure  show  the  mode  amplitude  time  histories  for  three 
wavelengths  in  the  system.  The  dashed  line  shows 
exponential  growth  at  the  unmodulated  growth  rate  of  the 
X.  = 100(j,  mode. 


AXk(cm) — - V(cm/sec) 


Dynamic  stabilization  of  long  wavelength  Rayleigh-Taylor 
modes  for  a 1 cycle  per  nanosecond  modulation.  The 
upper  curve  gives  the  average  shell  velocity  versus 
time  trajectory  for  the  runs  2B  and  2C.  3.5x10^  cm/sec 

was  achieved  without  instability  of  the  normal  R-T  type. 
The  three  wavelengths  in  the  lower  figure  all  lie  below 
the  dashed  line,  the  linear  growth  line  for  the 

= 300 p mode.  Note  that,  on  burnthrough,  the  outer 

surface  perturbation  is  still  only  0.1-0.2p. 
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